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We generate an explicit four-fold infinity of physically acceptable exact perfect fluid solutions 
of Einstein's equations by way of conformal transformations of physically unacceptable solutions 
(one way to view the use of isotropic coordinates). Special cases include the Schwarzschild interior 
solution and the Einstein static universe. The process we consider involves solving two equations 
of the Riccati type coupled by a single generating function rather than a specification of one of the 
two metric functions. 
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I. INTRODUCTION 

Perhaps the simplest of all procedures that one can 
think of for generating new exact solutions of Einstein's 
equations is the use of conformal transformations. Un- 
fortunately, when applied to vacuum, no new solutions 
emerge via this procedure [lj. However, when consider- 
ing static fluid spheres in isotropic coordinates, the seed 
metric is not vacuum, but an unphysical fluid with pres- 
sure but zero energy density and viable new solutions do 
indeed emerge. We now have a rather vast array of meth- 
ods for generating static fluid spheres 0], with many of 
the successful procedures relying on what amounts to the 
development of a linear equation of first order, for exam- 
ple and 0. In isotropic coordinates linear equations 
do not emerge directly @. Rather, specifying one of the 
two metric functions leads to a differential equation of 
the Riccati type. Here we do not specify either metric 
function but rather solve two equations of the Riccati 
type coupled by a single generating function. Whereas 
we are able to solve this system for a variety of generat- 
ing functions, we have found only one class of generating 
functions that gives rise to tractable and physically in- 
teresting solutions of Einstein's equations. 



II. GENERATING TECHNIQUE 

Consider spacetimes O that are static conformal trans- 
formations of seed metrics S where m 
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where dil 2 is the metric of a unit sphere (d9 2 + 
sin 2 {9)d(j) 2 ) and F is a freely specifiable function > 0. 
We suppose that the mathematical fluid associated with 
O is generated by streamlines of constant r, 9 and </>. 
Since this flow is shear free, the necessary and sufficient 
condition for ([T]) to represent a static perfect fluid is given 
by the Walker pressure isotropy condition [7[ 



where is the Einstein tensor. The energy density is 
defined by 8np(r) = ~G\ and the pressure by 8irp(r) = 
G r r and for we assume that p + p ^ [1] . Condition 
([3]), along with the definitions for p andp, is equivalent to 
the Tolman-Oppenhcimcr-Volkoff equation. The space- 
times S (for any constant F) do not represent physically 
acceptable static fluid spheres since they all have zero 
energy density. What we are interested in are spacetimes 
O that represent physically acceptable exact perfect fluid 
solutions of Einstein's equations. 



From condition wc find 

r$" + r$' 2 - $' + J(r) = 
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and = d/dr. The energy density is given by 
F" 3F' 2 2F' 

independent of and the pressure is given by 
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= fT r + 
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where $ and F are linked by |4]) . It is clear from ((6|) that 
F must have a local maximum at r — and so from ([5]) 
we must have J(0) = 0. 

Let us specify F and solve for $ from ([!]). The formal 
solution is given by 



b{r)dr + C 



where 



b 2 - 
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and C is a constant. Since equation © may be solved 
analytically only for certain J, we can ask what F gives 
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rise to this particular J? The answer follows from ([5]) 
and is given by 



where 



F = exp / b(r)dr + C 



Z' ~ b2 1 J n 
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2 r r 



(10) 



(11) 



and C is a constant. Alternatively, we can specify $ and 
solve for F. This is equivalent to (fTU)) with pip where we 
consider J generated from ([4|. Again, we can ask what 
4> generated a particular J. The solution is given by ([8|) 
with ([9]). In either case our ability to proceed revolves 
around our ability to solve Riccati equations of the type 
® and (HI]). 

In the usual way, equations (HJ) and (fTTj) can be written 
in linear form as 



■0 1 = 

r r 



and 



Jlp 
~2r 



(12) 
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where, up to a scale factor in t and a constant confor- 
mal factor, -0 = e* and -0 = 1/vF- Rather than specify 



F or $, here equations (|12|) and (|13|) are solved simul- 
taneously, coupled by the generating function J. The 
solutions for ip and ip are, of course, quite different in 
general and, as explained below, distinct even for J = 0. 



III. SOLUTIONS TO THE RICCATI SYSTEM 

Given a particular solution to a Riccati type equation, 
standard procedures [9j allow the construction of more 
general solutions. However, this procedure for generat- 
ing solutions usually starts from very simple particular 
solutions, and we have found no non-trivial known so- 
lutions applicable to the system (fT2")l and (fT3"j) . Rather, 
what we have done is to use the computer algebra system 
Maple [l(| to generate solutions to this system. 

To motivate our choice for J, first consider 

.4 



F = 



l + Br 2 Y 
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where A and B are constants and n is a ratio of integers. 
This gives 



2(2 - n)nB 2 r 3 
{1 + Br 2 ) 2 



J 



(15) 



which distinguishes two special cases for which J = 0: 
n = and n = 2. It is important to note that (fTi)) 



is but a special case that leads to (|15| . With J = 0, 
e* = C + Dr 2 where C and D are constants. The cases 
n = are physically unacceptable, since, as explained 
above, the associated energy density vanishes. All cases 
with n — 2 are conformally flat and so rep resent the well- 
known Schwarzschild interior solution 11]. A special case 
is given by C = 1 and D — B which is the Einstein static 
universe. (A cosmological constant A = AB/A can be 
introduced to give zero pressure.) 

Motivated by the foregoing, we have considered the 
generating functions 



J = 



2(2- n)nB 2 r 
(1 + Br 2 ) a 



(16) 



where a and b are integers, and have been able to solve the 
Riccati system (|12[) and (fT3)) analytically for the integers 
shown in Table [I] 

TABLE I: Analytic Solutions 



a 


b 


1 


1, 3, 5 


2 


1, 3, 5, 7 


3 


1, 3, 5 


4 


1, 3, 5,7 


5 


3, 5 


6 


3, 5, 7 



However, solving the differential equations does not 
mean that we can find a physically acceptable, or even 
tractable, solution to the Einstein equations. In some 
cases (e.g. a = 2, b = 5) we have been unable to construct 
the associated energy density and pressure simply due to 
the complexity of the background spacetime. In other 
cases (e.g. a = 1, b = 3) the most elementary physical 
requirements cannot be met (finite positive p and p at the 
origin r = with monotone decreasing values outward). 
Of the solutions represented in Table [J we have found 
only one case of physical interest: (a = 2, b = 3), that is, 



IV. PHYSICALLY ACCEPTABLE SOLUTIONS 

Now starting with (fT5|) from (fT2|) and (fl3|) we find 

e* = d(l + Br 2)(l+v'F)/2 + Ca(1 + Bj ,2 ){ l-VN)/2 

where 

A=2n 2 -4n+l, (18) 

and 



F = 



(C 3 (l + Br 2 ) n / 2 +C 4 (1 + Br 2 f- n l 2 f 



(19) 
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where the C x are constants. From (fT8|) we have n > 
1 + \/2/2 and n < 1 — \f2j2. Whereas the metric is 
remarkably simple, the resultant expressions for the en- 
ergy density and isotropic pressure are very long and not 
reproduced here. We resort to graphical demonstrations. 

The fact that the energy density is unaffected by $ is 
demonstrated in Figure [1] where we have varied C\ and 
Ci- Two sets of curves shown coincide. This degeneracy 
arises due to the fact that we have not set a scale for t. 
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have been generated by solving two equations of the Ric- 
cati type coupled by a single generating function rather 
than specifying one of the metric functions. Special cases 
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FIG. 2: Variation of C3 and C4. For all the curves C\ = 1, 
C2 = 1/4, n — 7/4 and B = 1. The nature of the curves can 
be recognized from Figure [1] For one set of curves C4 = 4 and 
C3 is varied: 1 for the thick curves, 3/4 for the regular curves 
and 1/2 for the dashed curves. For the other set C3 = 1 and 
C4 is varied: 4 for the thick curves, ~5.3 for the regular curves 
and 8 for the dashed curves. 



FIG. 1: Variation of C\ and C2. For all the curves C3 = 1, 
d = 4, n = 7/4 and B — 1. The top curve shows p/p(0). 
The curves which intersect the abscissa give p/p(0) and the 
horizontal curves show the square of the adiabatic sound 
speed. The configurations terminate and match onto vacuum 
at p = 0. Two sets of curves, which coincide, are shown. For 
one Ci = 1 and C2 is varied: 1/2 for the thick curves, 1/4 
for the regular curves and 1/8 for the dashed curves. For the 
other, C2 = 1/4 and Ci is varied: 1/2 for the thick curves, 1 
for the regular curves and 2 for the dashed curves. 



of these solutions include the Schwarzschild interior solu- 
tion and the Einstein static universe. The solutions are 
qualitatively similar to the Tolman IV solution (see for 
example [3| ) and so should be of interest for the study of 
internal properties of neutron stars [l2| . 



In Figure [5] we have varied C3 and C4. Again, two sets 
of curves shown coincide. This degeneracy arises due to 
the fact that the essential physics does not change under 
a constant conformal transformation. 



V. DISCUSSION 

An explicit four- fold infinity of new physically accept- 
able exact perfect fluid solutions of Einstein's equations 
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